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Abstract 

An analytical model for investigating vibra- 
tion and flutter of mistuned bladed disk assemblies 
is presented. This model accounts for elastic, 
inertial and aerodynamic coupling between bending 
and torsional motions of each individual blade, 
elastic and inertial couplings between the blades 
and the disk, and aerodynamic coupling among the 
blades. The disk is modeled as a circular plate 
with constant thickness and each blade is repre- 
sented by a twisted, slender, straight, nonuniform, 
elastic beam with a symmetric cross section. The 
elastic axis, inertia axis, and the tension axis 
are taken to be noncoincident and the structural 
warping of the section is explicitly considered. 

The blade aerodynamic loading in the subsonic and 
supersonic flow regimes is obtained from two- 
dimensional unsteady, cascade theories. All the 
possible standing wave modes of the disk and 
traveling wave modes of the blades are included. 

The equations of motion are derived by using the 
energy method in conjunction with the assumed mode 
shapes for the disk and the blades. Continuities 
of displacement and slope at the blade-disk junc- 
tion are maintained. The equations are solved to 
investigate the effects of blade-disk coupling and 
blade frequency mistuning on vibration and flutter. 
Results show that the flexibility of practical 
disks such as those used for current generation 
turbofans does not have a significant influence on 
either the tuned or mistuned flutter characteris- 
tics. However, the disk flexibility may have a 
strong influence on some of the system frequencies 
and on forced response. 
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Introduction 

Vibration, aeroelastic stability and response 
of bladed-disk assemblies employed in modern air- 
craft turbofan engines have been among the most 
difficult analysis problems encountered. The 
blades are thin, and can have a considerable amount 
of built-in pretwist. As a consequence of pretwist 
and rotation, the blade natural modes involve coup— 
ling between the in-plane (the plane of rotation) 
and out-of-plane bending and the torsional motions. 
Another consequence of .rotation is. that the blades 
experience steady-state equilibrium loads which in 
turn require either an explicit or implicit con- 
sideration of geometric nonlinear theory of elas- 
ticity to adequately model the blade behavior. 

Additionally, the individual blades are struc- 
turally coupled through the disk and often through 
a variety of other connecting parts such as 
shrouds, dampers, and lacing wires. These problems 
are further compounded by the presence of small 
differences between the individual blades, known 
as structural mistuning. 
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Unlike the case of a fixed-wing, the flow Mach 
number varies along the blade span and involves 
unsteady subsonic, transonic, and supersonic flow 
regimes. The flow field is three-dimensional and 
also involves aerodynamic interference between 
adjacent blades within a stage as well as from 
upstream and downstream stators and rotors. 
Furthermore, there may be small differences in the 
flow between the individual blades, known as aero- 
dynamic mistuning. 

Because of the complexities mentioned above, 
the research work in the subject area has been 
proceeding with simple but progressively refined 
mathematical models. Vibration of the bladed-disk 
assemblies with structural coupling have been 
studied in many publications, for example Refs. 1 
to 4. Flutter and response using single degree-of- 
freedom blade models without structural coupling 
between blades, including mistuning, are addressed 
in Refs. 5 to 8. Refs. 9 to 14 have considered 
flutter of tuned assemblies with rigid disks. A 
finite element formulation was presented in Ref. 

15 to study flutter of bladed-disk systems. 

Refs. 16 and 17 investigated the effects of mis- 
tuning on cascade flutter and response in incom- 
pressible, subsonic and supersonic flow regimes 
with a rigid disk and a typical section model in 
which each blade has two degrees-of-freedom, one 
bending and one torsion. Effects of blade mis- 
tuning on flutter were also investigated in Refs. 

18 to 19. Furthermore, the single degree-of- 
freedom typical section model used in Ref. 5 was 
modified in Ref. 20 to include disk flexibility 
and a sub- sonic aerodynamic theory was used to 
investigate flutter. The two degree-of-freedom 
typical sec- tion model used in Refs. 16 and 17 
was refined in Ref. 21 to include disk flexibility 
and to study flutter and response. 

.The typical section model used in Refs. 16, 

17 and 21 was found to be adequate to elicit phys- 
ical understanding of mistuning effects. Paramet- 
ric studies indicated that structural mistuning 
has the potential to significantly raise the flut- 
ter speed of an advanced turbofan. This potential 
may be utilized to eliminate the shrouds which 
have an adverse effect on aerodynamic performance. 
To obtain more accurate flutter boundaries, includ- 
ing the effects of mistuning, the typical section 
model was replaced by a beam model with a rigid 
disk in Ref. 22. The beam model was further ex- 
tended in Ref. 23 to include blade sweep in order 
to investigate the possibility of classical flutter 
in advanced turboprops. 

A logical extension to the present state of 
the literature is to refine the structural model 
used in Ref. 22 to include disk flexibility by 
treating the disk as a circular plate. The purpose 
of the research in the present paper is to combine 
such a structural model for the disk with a beam 
model for each blade in order to study the effects 
of coupling between disk and blade modes and of 
blade frequency mistuning on vibration and flutter.. 

The beam model used in Ref. 22 is similar to the 
one used in Ref. 24, but there the effect of warp- 
ing was only partially considered in an implicit 
manner. In Ref. 25, the present authors found 
that the partial inclusion of warping in calculat- 
ing uncoupled torsional frequencies and mode shapes 
of twisted rotating turbofan blades is not adequate 
and, hence, they modified the model to include all 


the warping effects. This modified beam model will 
be used herein. The unsteady, two-dimensional, 
cascade, aerodynamic loads are calculated by 
Smith's theory'*®' in subsonic flow and Adamczyk 
and Goldstein's theory^*') in supersonic flow 
with subsonic leading edge; 

The equations of motion are derived by using 
an energy method in conjunction with assumed mode 
shapes for the disk and for each blade. Continuity 
of displacements and slopes at the blade-disk junc- 
tion is maintained. The resulting equations are' 
cast in the form of a standard complex eigenvalue, 
problem from which the vibration and aeroelastic 
stability characteristics are determined. A com- 
puter program is developed to form and to solve 
the complex eigenvalue program. In order to have 
checks at least for some special cases, the 
equations and the computer program developed in 
Ref. 22 form special cases of the present ones. 

This paper presents an analytical model, the 
development of the equations of motion and computer 
program, and some selected results for an advanced 
fan stage. 

Theory 

The motion of a bladed-disk assembly with or 
without mistuning can be expressed in terms of 
either a standing or traveling wave form. Since 
the stability is independent of the wave represen- 
tation, consideration of either of the representa- 
tions is adequate to describe the motion. In most 
of the published literaturel®’ 1 ®’* 1 ’**) dealing 
with blade flutter of turbomachines, the traveling 
wave modes (also called interblade phase angle 
modes) are preferred because the unsteady aero- 
dynamic loads are expressed in these modes. This 
is in contrast to the conventional flutter analysis 
of nonrotating wings and vibration of disks in 
which the standing wave modes are used. However, 
the present formulation uses both types of modes 
because the development includes both the blade 
and the disk flexibility and also is based on the 
development in Ref. 22. Hence, it will be neces- 
sary to understand both the traveling and standing 
wave mode representations and the relation between 
them. 

To understand the traveling wave modes, con- 
sider the vibratory bending motion of each blade 
of an N-bladed rotor with a rigid disk; For illus- 
tration purposes, consider only one bending mode 
for each blade. Then the cascade has N degrees- of 
freedom. The motion of the system can be expressed' 
in terms of traveling wave modes each of which is 
characterized by the phase angle B r between 
the adjacent blades. The angle is commonly called 
the interblade phase angle. The phase angle is 
restricted to N discrete values, 8 r = 2nr/N where 
the wave number index r = 0, 1, 2, ... { N— 1 ) . In 
each of these modes, all blades move with the same 
amplitude and every adjacent pair moves with the 
same interblade phase angle. 

To understand the standing Wave modes and the 
relation between the traveling and standing wave 
modes, consider the disk bending deflection perpen- 
dicular to the plane of the disk in the traveling 
wave form as 
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Motion of a Tuned Bladed-Disk System 
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R H 


_■ Kr* 

Zw R(l 'd ) u dar e 
r=0 



(1) 


Theoretically, the wave numbers for an elastic disk 
can have values up to infinity. However, if N 
blades are attached to the disk and if the number 
of traveling wave modes are restricted to N, then 
the wave number index is also limited to (N-l). 
Thus, the disk with the N attached blades can be 
characterized by N traveling wave modes. Assum- 
ing simple harmonic motion, the real part of 
Eq. (1), can be rearranged to give the disk deflec- 
tion in the standing wave form as 


u d 


R ( n d) t u dacO + u dacN/2 


cos 


N 

2 ♦ 


Following common practice in the literature, 
the disk motion is expressed in standing wave modes 
and the blade motion is expressed, as in Ref. 22, 
in traveling wave modes. For the disk, using a 
single quadratic deflection shape along the radius 
for all the possible nodal diameter modes, the 
expression for the motion is given by Eq. (2) with 


R(n d ) = n d (3) 

For the blades, assume the number of modes retained 
in the plane of rotation, in the plane perpendicu- 
lar to the plane of rotation, and in torsion are 
Mh, Mg, and M a , respectively. Then, the 
motion of the sth blade can be represented as 


+ 



< u dacn cos n * 



+ u dasn sin n e ° ( 2 ) 

where N D = (N-l)/2 for odd N and N 0 = ( N-2 ) / 2 
for even N, and n is the nodal diameter index. 
Thus, the disk with the N attached blades can be 
characterized by 0,1,2,... Nq and N/2 nodal 
diameters. The N/2 nodal diameter (second term in 
the square bracket of Eq. (2)) is present only for 
even N. Comparing Eqs. (1) and (2), it is clear 
that the traveling wave mode with the zero inter- 
blade phase angle, ^ da Q, corresponds to the zero 

nodal diameter mode, q, that the traveling wave 

mode (with the interblade phase angle = ir),- %aN/2’ 
corresponds to the N/2 nodal diameter mode, and the 
traveling wave modes u dflr and u d a ( N— r ) com ^ ne to 

form two standing wave modes, ^ dacr and u^ . 

Thus, either a traveling or a standing wave mode 
approach can be used to describe the motion of a 
disk and the same is true of the blade or of the 
bladed-disk system (see also the discussion in 
Refs. 14, 16, 22 and 28). The above relationship 
will be used in deriving the blade aerodynamic 
loads due to the disk motion. 

Coordinate Systems 


Coordinate systems employed in deriving the 
equations of motion of the sth blade are shown in 
Figs. 1 and 2. The axis system Xo, Yn, Z n 
rotates with a constant angular velocity n about 
the Xjj-axis. The Yjj-axis coincides with the 
undeformed elastic axis of the blade. The prin- 
cipal axis, x, of the blade cross section at any 
point on the elastic axis is inclined to the 
X n -axis by an angle s as shown in Fig. 2. The 
disk elastic deformation, u d , and the blade defor- 
mations, u, v, w, and a, translate and rotate the 
xyz system to the x 3 y 3 Z 3 system. 
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In Eqs. (4) to (6), the complex constants h^, g 

and a - which are commonly used in fixed wing and u 
rotary s wing aeroelasticity are expressed in terms 2 



of the new constants, h ar j, g ar j, and a ar j and 

interblade phase angle since the aerodynamic loads 
are expressed in terms of the new constants. The 
total number of degrees of freedom for the bladed 
disk system is N times (1 + M^ + M^ + .M a ). The 

function f (g) in Eq. (5) is determined by enforc- 
ing the continuity of displacement and slope at the 
blade-disk junction point 
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The nonrotating frequency, for each mode is 
obtained from 

2Fg + (2Fg + F 2 ) cos kg cosh k^ . -• 

+ F^ ypg" sin k^ sin hk^ = 0 


( 13 ) 


(14) 


In calculating the mode shapes and frequencies, 
average values of actual blade shape are used for 
L/c, c/t m and £ in Eq. (13). 

Motion of Mistuned Bladed-Disk System 


These two conditions are satisfied if 

f c<") - sj; <1 4 2 " <9> 

In Eqs. (4) and (5), the standard nonrotating 
orthogonal normal modes in bending for a beam with 
fixed-free boundary conditions are given by 


Wj.(n) = cosh (PjTT) - cos (p^n ) 


(cos Pj + cosh Pj) 
(sin pj + sinh pj) 


[sinh( Pj . n ) - sin(pjn)] 


( 10 ) 


where Pj is calculated from 

cos Pj cosh pj + 1 = 0 (11) 

In Eq. (6) the orthogonal nonrotating modes in 
torsion for a uniform twisted beam with fixed-free 
boundary conditions are obtained from Ref. 25 and 
are 

— — ^1 — 

Aj(u) = [(sinh kjn - y- sin kgn) 


p 

(kj sinh kj + k^kg sin kg) 

2 2 " 

(kj cosh kj + kg cos kg) 


(cosh k^n - cos kgn) ] 
( 12 ) 


where 


In the present formulation we consider mis- 
tuning only due to nonidentical structural proper- 
ties of the blades, since the disk thickness is 
assumed to be uniform. Because of the mistuning 
there will be coupling between various traveling 
wave modes or equi valently between the various 
standing wave modes. However * the general motion 
of an arbitrarily mistuned, system can be expressed 
as a linear combination of the motions in either 
traveling or standing wave modes. The disk motion 
in Eq. (2) is already expressed in terms of the 
possible nodal diameter modes. The motion of the 
sth blade is obtained by simply superposing the 
blade deflections (in Eqs. (4) to (6)) in the rth 
traveling mode from r = 0 to r = N-l. 

Aerodynamic Model 

For a rigid disk, the motion dependent aero- 
dynamic lift and moment per unit blade span are 
expressed in Ref. 22 in terms of four coefficients 
*hhr» ®-hor> ^ahr> f aa r an< l blade motions 
which are written in a traveling form. These co- 
efficients were calculated by using Smith's theory 
in subsonic flow and Adamczyk and Goldstein's 
theory in supersonic flow with a subsonic leading 
edge. At any radial station, the relative Mach 
number is a function of in-flow conditions and the 
rotor speed. Most current fan designs have super- 
sonic flow at the blade tip and subsonic flow at 
the blade root. Consequently, some region of the 
blade encounters transonic flow. Since the above 
theories are not valid in the transonic region, 
the subsonic theory with M ?f f = 0.9 for stations 
in the range 0.9 < M f f < I and the supersonic 
theory with M e jf = 1.1 for stations in the range 
1.0 £ M e ff £ 1.1 were used. 



For a flexible disk, the lift and moment ex- 
pressions can be divided conveniently into two 
parts: one due to blade motion and the other due 

to disk motion. The first part directly flows 
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from Ref. 22. The second part can be obtained in 
two steps. First, express the lift and moment due 
to disk motion in the traveling wave form repre- 
sented by Eq. (1) in terms of the four coeffi- 
cients, mentioned earlier. Second, by taking the 
real part of the resulting expressions, by re- 
arranging terms, and by assuming simple harmonic 
motion, derive the expressions for lift and moment 
per unit span due to disk motion represented by 
Eqs. (2) and (3) (see also refs. 14 and 28). Com- 
bining these two parts, the lift and moment expres- 
sions per unit span of the sth blade can be derived 
and are 
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where 


£ hhln ” 7 |*hhn + £ hh(N-n)j 
1 hh2n = " 7 [ £ hhn " £ hh(N-n)j 


(17) 


<*ha2n u dasn + £ haln u dacn> sin 


1 T 


(15) 


£ aa2n 7 aan £ aa(N-n)J 

The coefficients £ hhr , t h r , i hr , and i aar 
are calculated for each vaTue of r for specified 
values of Mach number, reduced frequency, s/c, s, 
and a. Then, the corresponding coefficients 
*hhln, *hh2n»---> *aa2n for each value of the 
nodal diameter index n are obtained from 
Eq. (17). Note that the subscripts r and n 
represent the interblade phase angle and nodal 
diameter number, respect? vely. 


6 



Structural Model 


and are 


The bladed-disk system is idealized as a uni- 
form circular disk with constant thickness to which 
the blades are rigidly attached. Only the disk out 
of plane bending motion is considered. The blade 
model is the same as that used in Ref. 22 except 
for the fact that the warping of the blade cross 
sections is explicitly considered herein. Conti- 
nuity and slope at the blade-disk junction is en- 
forced as described in the previous section. 

Equations of Motion 

The strain energy expression for a circular 
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where 


, D = E d td /12 


(19) 


The stress components o r and are obtained 
by solving the following equations 
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( 21 ) 


The strain energy of the s tfl blade follows from 
Ref. 22 and is 
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The explicit consideration of warping introduces 
two additional coefficients Ci and C 2 and 
modifies the expression for J: 


Cj = x 2 dx dz 
C? = yy x z dx dz 

■//to - i> 2 * (« * Hn 1” i 


s. (23) 


The other coefficients in Eq. (22) are defined in 
Ref. 22. 

The kinetic energy of the disk is 
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The kinetic energy of the s^ blade also follows 
from Ref. 22 and is. • 1 
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The Coriolis and the inertial warping terms are 
considered to be small and are neglected. 


The total strain and kinetic energies of the 
system are 
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The generalized aerodynamic forces are calculated 
from the following virtual work expression 
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By substituting Eqs. (26) and (27) into Lagrange's 
equation, the equations of motion of the system 
are obtained. The final equations are nondimen- 
sionalized and cast in a standard eigenvalue form 
as 
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[Q] = [M] + [ED] [A] [ED ] -1 


(X) = [ED] {Y} 


(X) 


[ED] 


<v 


<v 

<v 


l< Y o » / 

< x l> 

; { Yi = 

K > 

f * 

\ {X N-1> 


{Y (N-1) } ] 


[ 1 ] [ 0 ] 

[0] [E] 


(28) 


(29) 


<v 


u dacO ' 


*s. \ j 

' W arl 

u dacN 

2 


*s 2 


*ar 2 

u dacl 




i ; 

: 

> 

1 ? sl 


j 9 ar l 

“dacNjjl 

; {X $ i = 

Us 2 

; cy = 

j ^ar 2 

u dasl 1 


“si 1 


' “arl 

• 


“s 2 


“ar 2 

“dasN Dj 




l'- 


The matrices [P], [M] and [A] are the stiffness, 
mass and aerodynamic matrices of the system. The 
matrix [E] is defined in Ref. 22. 


Results and Discussion 


Solution 


The flutter boundaries are obtained by solving 
the standard complex eigenvalue problem, Eq. (28). 
The relation between the frequency w 'and y is 



i v 


(30) 


Flutter occurs when v > 0. Vibration characteris- 
tics are obtained by setting the aerodynamic matrix 
[A] to zero. 

A computer program was written to assemble and 
solve the generalized eigenvalue problem. As 
mentioned earlier, the formulation and the program 
have very wide scope for investigating the effects 
of disk flexibility and blade frequency mistuning 
on both vibration and flutter of bladed disk assem- 
blies. Because of space limitations, only selected 
results are presented. 


Vibration 


An advanced unshrouded fan stage (aspect 
ratio = 3.3) representative of a next-generation 
fan was chosen for analysis. This fan was analyzed 
in Ref. 22 by treating the disk as rigid. The disk 
of this fan consists of a thick hub, a thinner web 
and a short thick rim to house the blades. Accord- 
ingly, a disk model which consists of three concen- 
tric connected nonuniform annuli would be required 
for accurate modeling. However, to reduce the 
complexity of the model, the disk of the fan is 
idealized as a circular disk with uniform thick- 
ness. To investigate the effects of disk flexi- 
bility on both vibration and flutter, four values 
for disk thickness are considered. For descriptive 
purposes, the disk with thicknesses 1.891, 0.1891, 
0.0946 and 0.0473 m are designated as very rigid, 
rigid, flexible, and very flexible disks, respec- 
tively. The equivalent thickness for a real prac- 
tical fan may lie between the values for the rigid 
and flexible disk. The other properties of the 
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disk are'listed in Table 1. The blade properties 
are the same as those listed in Ref. 22 and hence 
they are not .repeated here. The design speed of 
this 28-bladed fan stage is 4267- rev/min. 

The vibration analyses are performed by using 
six modes for each of the. 28 blades and all pos-. 
sible nodal diameter modes for the disk. The blade 
modes consist of two modes each in the plane of 
rotation, in the plane perpendicular ,'to the plane 
of rotation, and in torsion. The total number of 
degrees of freedom is 28 x 6 + 28 = 196; 

The blade frequencies with disk thickness 
tq = 1.891 m are virtually the same as the blade 
cantilever frequencies. Thus, the disk with 
tq =1.891 m behaves virtual ly like an infinitely 
rigid disk. Since the very rigid disk case was 
analyzed in Ref. 22 for vibration and flutter, the 
results are not repeated herein. The variation of 
the frequency ratio with nodal diameter index is 
i 1 lustrated. in Figs. 3 to 5 for the other three 
values of. the disk thickness. . The frequency ratio 
is obtained by dividing the frequency by an arbi- 
trarily chosen reference value, u 0 = 200 Hz. In 
these. figures the frequency ratios are designated 
as F,| C,vT and D indicating, that the blade flat- 
wise, chordwise, torsion and disk motions, respec- 
tively,: are predominant "in the eigenvectors. The 
additional numbers 1 and 2 in the designations 
denote first and second modes. The blade canti- 
lever, f requencies are shown by dark symbols. Based 
on these frequency ratios and the corresponding 
eigenvectors, which are not. shown because of space 
limitations, several interesting observations are 
made. (1) For a 28-bladed system, the identifica- 
tion of the. predominant motion in the eigenvector • 
was found to be a tedious, job. This problem was 
solved by analyzing a system with 4, 8, 12 and 16 
blades since the identification of predominant 
motion is easier for fewer blades. (2) For each 
nodal diameter, except for zero and fourteen, each 
eigenvalue represents a pair of. modes and each 
mode in the. pair is a linear combination of the 
sin n^ s and cos ni|) S . ’ Furthermore, the 
mode's Tn each pair are orthogonal as expected for 
a tuned -system. The eigenvalues for the zero and 
the 14th nodal diameter represent only one mode. 

In the zeroth nodal diameter .mode all the blades 
move in phase and in the 14th nodal diameter mode 
the phase difference between one blade and its 
neighbor is 180°.. (3) Because of the pretwistof 
the blade, there is axoupling between blade flat- 
wise and chordwise motions and between blade chord- 
wise and disk motions. For the disk with 
tq = 0.1891 m, it can be seen from Fig. 3, that 
the coupling between blade chordwise and the disk 
motions is stronger than that between blade flat- 
wise and disk motions'. This is because the disk 
frequency is- in the vicinity of ' the blade second 
chordwise frequency. This is contrary to assump- 
tions often made in the published literature deal- • 
ing with bladed-disk vibrations. (4) Comparing 
the results in Figs. 3 to 5, the disk frequency 
decreases with a decrease in disk stiffness, as 
should be expected. Furthermore, the coupling 
between the disk and blade motions and blade to ; 
blade motions increases as the disk becomes thin- 
ner.' Alternatively, one could say that'more of 
the disk modes couple with the blades modes as the 
disk flexibility decreases:' : (5), For the disk with 
tp = 0.1891 m, it can be seen from Fig. 3 that 
tne coupling between the disk and blades for nodal 
diameters higher than four is almost negligible. 


The effect of alternate mistuning on vibration 
is investigated. The method used here to vary the 
frequency from bl ade-to-blade is-simply to vary 
either the blade torsional stiffness GJ or the 
bending stiffness Ef^lxx • For alternate mistuning 
the torsional stiffness of the odd numbered blades 
is increased by 10 percent and that of the even 
numbered, blades is decreased by 10 percent from 
that, of a tuned blade system.' The result is a 
frequency variation Xf approximately 7 percent. 

When this mistuning is introduced, the; T-l and T-2 
modes in Figs. 3 to 5 are affected. : The frequency 
ratios for the T-l mode' alternates 'between 1.06 
and 1.13 and that for the T-2 mode between 2.75 
and 2. 93. ; The effect of mistuning in the other 
frequencies is- very small.. 

Flutter of a. Tuned Bladed-Disk -. 

The advanced fan stage analyzed for vibrations 
is again considered to investigate flutter behavior 
including disk flexibility. The solution technique 
used herein to calculate the flutter boundary is 
the same as that described in Ref. 22. The cor- 
rectness, of the program was checked by comparing, 
the flutter speed for the case with tq = 1.891 m ,. 
with that for the rigid disk case in Ref. 2,2 in- 
cluding all the warping effects. Next, the" effect 
of disk flexibility on flutter speed is illustrated 
in Fig. 6. The fan stage' is stable, if the operat r 
ing point is to the left of the boundary. For a 
! disk with tq = 0.1891 m the flutter boundary 
1 coincides with that for the tuned rigid disk 
( tq = 1.891 m), indicating that this amount of 
disk’ flexibility does not have any influence on 
flutter. For a disk thickness tq =0.0946 m, . 
the flutter speed is very slightly reduced. This 
reduction in flutter speed is more significant, 
when the disk thickness is reduced to 0.0473 m. 

This degradation in flutter speed is expected 
because the coupling between disk. and' blade becomes 
very strong and reduces blade frequencies (see 
fig. 5). Additional ly, the flutter mode is changed' 
from a predominantly torsional mode for the rigid 
disk to a predominantly bending mode for a very 
flexible disk. It should be mentioned. that there 
is ho blade bending-torsion coupling inherent in 
the blades themselves because the center of gravity 
and elastic axis of any blade are coincident. 
However, coupling between bending and torsional 
motions does, occur because of the disk flexibility. 

Furthermore, there is aerodynamic coupling between 
bending and torsional motions. Thus, flutter.' in- 
volving a. predominantly bending mode is not sur- 
prising (see also, Ref. 16). 

Flutter of a Mistuned Bladed^Disk 

Previous publications^, 17, 22) have ^ 0 ^ 
that the torsional . frequency mistuning. could have a 
significant stabilizing effect on flutter. The 
effect of mistuning is further investigated herein 
with theflexible disk model and some selected 
results are presented in Fig. 6. Alternate blade 
torsional frequency mistuning of approximately 7 
percent has significantly increased the flutter 
speed, for .the cases with tq = 0.1891, and 0.0946 m. 
This finding is in agreement with that in Ref. 22. 
Again, .this amounts to saying that the disk flexi- 
bility does not have much influence on the flutter 
speed even in the presence of mistuning. Contrary 
to these two cases, there is no appreciable im- 
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provement on flutter speed due to torsional fre- 
quency mistuning for the very flexible disk. This 
might be expected because the predominant motion 
in the flutter mode is bending rather than torsion. 
Next it is natural to ask whether an alternate 
bending frequency mistuning has any influence on 
flutter speed. To investigate this, the bending 
stiffness EI m of the odd numbered blades was 
increased by 10 percent and that of the even num- 
bered blades was decreased by 10 percent from that 
of a tuned blade. The result is a frequency vari- 
ation of approximately 2.6 percent. With this 
mistuning the flutter boundaries for the very 
flexible disk is shown again in Fig. 6. Evidently, 
there is a small increase in flutter speed. How- 
ever, this increase in flutter speed is not as 
great as that due to alternate torsional mistuning 
for the rigid disk. Again, it is concluded that 
mistuning has potential to increase the flutter 
speed even in the presence of mechanical coupling 
between blades introduced by disk flexibility. 
Qualitatively, these findings are in agreement 
with the corresponding ones in Ref. 21 where a 
modified typical section model was used. However, 
it should be emphasized that the flexibility of 
the disks which are employed for the current gen- 
eration of turbofans does not have a significant 
effect on flutter speed. Still, the flexibility 
effects should be included in calculating forced 
response which is sensitive to small changes in 
system frequencies. 

Concluding Remarks 

An analytical model and an associated computer 
program to investigate the effects of disk flexi- 
bility on vibration and flutter characteristics of 
both tuned and mistuned bladed-disk assemblies are 
developed. From the parametric studies the follow- 
ing conclusions are reached. 

(1) The flexibility of the disks which are 
employed for the current generation turbofans does 
not have much influence on either tuned or mistuned 
bladed-disk system flutter characteristics. Still, 
its effect may be significant on forced response 
which is sensitive to small changes in system 
frequencies. 

(2) The flexibility of a very thin disk re- 
duces flutter speed and changes the predominant 
character of the flutter mode from torsion to 
bending. 

(3) Because of the blade pretwist, chordwise 
motion of the blades is more strongly coupled to 
the disk motion than is the blade flatwise motion. 
This is contrary to assumptions often made i n the 
published literature dealing with bladed-disk 
vibrations. 

(4) A moderate amount of blade frequency 
mistuning has significant potential to increase 
the flutter speed, even in the presence of the 
mechanical coupling between blades that is intro- 
duced by disk flexibility. 
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TABLE I. - DISK PROPERTIES 




Figure 1. - Blade coordinate systems before 
and after deformation. . 
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